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Abstract 

Based on a quantum Langevin equation and its corresponding Hamiltonian within a c-number 
formalism we calculate the vibrational dephasing rate of a cubic oscillator. It is shown that leading 
order quantum correction due to anharmonicity of the potential makes a significant contribution 
to the rate and the frequency shift. We compare our theoretical estimates with those obtained 
from experiments for small diatomics N2, O2 and CO. 
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I. INTRODUCTION 



A molecule in a liquid undergoes random frequency fluctuations due to the stochastic 
forces imparted by the surrounding medium. The correlation between this frequency fluc- 
tuations results in vibrational dephasing. The problem has received wide attention both 
from theoretical and experimental point of view over the last couple of decades. Several 
approaches to understand the basic nature of vibrational dephasing have been made. These 
include notably binary collision theorjJ- in which fluctuations are taken into account in 
terms of collisional events, hydrodynamic mode Pi^i'^i^ relating random force and shear vis- 
cosity of the liquid, Hamiltonian theory^ in terms of normal mode analysis and numerical 
simulation3^i2i2ii2iii using various molecular potentials. A key element of these formulations 
is the realization that vibrational dephasing owes its origin to cubic anharmonicity of the 
vibrational potential. In the present paper we attempt to explore further this issue within 
a quantum mechanical content. 

A good number of approaches to vibrational dephasing make use of generalized Langevin 
equation that governs the dynamics of the system by an inflnite number of harmonic oscilla- 
tors coupled linearly to the system. Very recently based on a coherent state representation 
of noise operator and a positive deflnite Wigner canonical thermal distributioni^ of bath os- 
cillators a c-number quantum Langevin equatio n^°i^^i^^i^^i^^i^^i^^ in the context of rate theory 
and stochastic processes has been proposed. An important offshoot of this formulation is 
that it takes care of quantum correction due to nonlinearity of the system potential order 
by order. It would thus seem that one should be able to analyze the vibrational dephasing 
rate quantum mechanically for arbitrary noise correlation and temperature and explore the 
role of this quantum correction to anharmonicity. This is the main purpose of this paper. 
In what follows we construct a c-number Hamiltonian within a normal mode description 
and estimate the dephasing rate by calculating the effective correlation between the relevant 
c-number co-ordinates of anharmonic origin. Since for a small diatomic molecule like N2, 
the fundamental vibrational frequency ujq is so high (2326 cm~^) that the ratio tiWo/kT is 
as large as 43.52 even at, say, 77 ^K, it is imperative that quantum nature of the oscillator 
molecule is signiflcant in estimating the dephasing rate in the harmonic as well as in the 
anharmonic regime. With this end in view we examine the vibrational dephasing rate to 
estimate the anharmonic quantum correction to this rate and its variation with temperature 
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away from critical or triple point for three widely studied diatomics A^2, O2 and CO to allow 
ourselves a fair comparison with experiment o^^i^^i^^i^^i^^i^^i^^ . 

The outlay of the paper is as follows: In Sec. II we introduce the quantum Langevin equa- 
tion and its Hamiltonian counterpart within a c-number normal mode description. Since the 
cubic nonlinearity gives rise to a leading order contribution to dephasing rate we estimate 
the quantum vibrational dephasing rate for a cubic oscillator in Sec. III. The quantum cor- 
rections due to nonlinearity of the system potential is calculated explicitly in Sec. IV. Sec.V 
is devoted to the results obtained theoretically for three diatomics N2, O2 and CO which 
are compared with experiments. The paper is concluded in Sec. VI. 

II. C-NUMBER QUANTUM LANGEVIN EQUATION AND NORMAL MODE 
TRANSFORMATION 

We consider a particle of mass /i coupled to a medium comprised of a set of harmonic 
oscillators with frequency Ui. This is described by the following Hamiltoniaii^: 



Here q and p are co-ordinate and momentum operators of the particle and the set {xi,pi} 
is the set of co-ordinate and momentum operators for the reservoir oscillators of mass 
coupled linearly to the system through their coupling coefficients q. The potential V{q) is 
due to the external force field for the Brownian particle. The co-ordinate and momentum 
operators follow the usual commutation relations [q,p] = ih and [xi.Pj] = iMij. 

Eliminating the reservoir degrees of freedom in the usual way we obtain the operator 
Langevin equation for the particle, 



H 






(2.2) 



where the noise operator F[t) and the memory kernel 7(t) are given by 



j 




%(0) + g(0) 




■2 COSUJjt + 



■Pj{0) sin Ujt 



(2.3) 
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and 



N ^2 



7it) = y2^^cosujjt , (2.4) 



,=1 ^^-^J- 



The Eq. fl2.2|) is the well known exact quantized operator Langevin equation for which the 
noise properties of F{t) can be derived by using a suitable initial canonical distribution of 
the bath co-ordinate and momentum operators at t = as follows; 



{m)QS = (2.5) 
l{{mm)QS + {F{t')F{t))Qs} = ^X^^/k^Jcoth^'jcosc^.^t- 0(2.6) 

where {..■)qs refers to quantum statistical average on bath degrees of freedom and is 
defined as 

{d)Qs = T^Q^^P(-^^^th/A:^T) ^2.7) 
TTexp{-Hbi,th/kBT) 



for any operator d{{{mjijj /cj)xj + q},{pj}) where -f^bath(Z]iIi(P?/2"^i + mi/2{ujiXi + 
;j^^g)^)) at t = 0. By Trace we mean the usual quantum statistical average. Eq. (j2.6|) is 
the fluctuation-dissipation relation with the noise operators ordered appropriately in the 
quantum mechanical sense. 

To construct a c- number Langevin equatio n^°i^^i^^i^^i^^i^^i^^ we proceed from Eq. ()2.2|) . We 
first carry out the quantum mechanical average of Eq. ()2.2|) 

^^m) + ^^ fdt'^it - t'mt')) + {V'm = {Pit)) (2.8) 
Jo 

where the quantum mechanical average (...) is taken over the initial product separable 
quantum states of the particle and the bath oscillators at t = 0, |0){|ai)|«2) • • • I^Af)}- 
Here 10) denotes any arbitrary initial state of the particle and \ai) corresponds to 
the initial coherent state of the i-th bath oscillator. \ai) is given by !«,) = 
exp(— |a;jp/2) ^^^Q(a"YA/nJ)|nj), being expressed in terms of the mean values of the 
shifted co-ordinate and momentum of the i-th oscillator, {{niiuf /ci){xi{0)) + {q{0))} = 



^^h/2miUJi{ai + a*) and (pj(0)) = i^hmiUJi/2[al — a^), respectively. It is important to note 
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that {F{t)) of Eq. (j2.8j) is a classical-like noise term which, in general, is a non-zero number 
because of the quantum mechanical averaging and is given by {{F{t)) = f{t)); 



m 



E 



:%(0)) + (g(0)) 



cos Ujt + 



rrijuj] 



rrijUJj 



{pj{0)) sin Ujt 



(2.9) 



It is convenient to rewrite the c-number equation ()2.8|) as follows; 



mt))+f^ / dt'^it-t'){m) + {v'm = fit) 



(2.10) 

To realize /(t) as an effective c-number noise we now introduce the ansat a^^i^^i^'^i^^i^'^i^^i^^i^^ 
that the momenta {pj{0)) and the shifted co-ordinates {(mjc<j|/cj)(xj(0)) -|- (g(0))}, {pj} of 
the bath oscillators are distributed according to a canonical distribution of Gaussian form 
as 



J J 



Vj =Afexp 



(2.11) 



so that for any function of the quantum mechanical mean values of the bath operator 
Oj{{pj{0)), {{nijUJj /cj){xj{0)) + (g(0)))) the statistical average (. . .)s is 



{0,)s = I O, V, d{p,{0)) d{imju]/c,){x,m + (g(0))} . (2.12) 

Here Uj indicates the average thermal photon number of the j-th oscillator at temperature 
T and nj{oJj) = l/[exp {tvjjj/ksT) — 1] and M is the normalization constant. 

The distribution ()2.11|) and the definition of statistical average ()2.12j) imply that f(t) 
must satisfy 



{fit))s = (2.13) 

and 

{mm)s = IY1 fcoth ^) cos.;,(t - t') (2.14) 

That is, c-number noise f{t) is such that it is zero-centered and satisfies the standard 
fluctuation-dissipation relation (FDR) as expressed in Eq. fj2.6|) . It is important to emphasize 
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that the ansatz ()2.1ip is a canonical thermal Wigner distribution^^ for a shifted harmonic 
oscillator which remains always a positive definite function. A special advantage of using 
this distribution is that it remains valid as pure state non-singular distribution function at 
T = 0. Furthermore, this procedure allows us to bypass the operator ordering prescription 
of Eq. (j2.6|) for deriving the noise properties of the bath in terms of fluctuation-dissipation 
relation and to identify f{t) as a classical looking noise with quantum mechanical content. 
The procedure has been used by us in several recent context o^^i^-^i^^i^^i^^i^^i^^ . 

We now return to Eq. ()2.10p to add the force term V'{{q)) on both sides of Eq. ()2.10p and 
rearrange it to obtain 



fiq = p (2.15) 

P = - f dt'^it - t')p{t') - V\q) + /(t) + Q{t) (2.16) 
Jo 

where we put {q{t)) = q{t) and {p{t) = p{t) for notational convenience and 

Q{t)=V'{m-{V\q)) (2.17) 

represents the quantum correction due to the system degrees of freedom. Eq. ()2.16|) offers 
a simple interpretation. This implies that the quantum Langevin equation is governed by a 
c- number quantum noise /(t) originating from the heat bath characterized by the properties 
(j2.13|) and (j2.14j) and a quantum fluctuation term Q{t) characteristic of the non-linearity of 
the potential. 

Referring to the quantum nature of the system in the Heisenberg picture, one may write. 



q(t) = q + 5q (2.18) 
p{t) = p + 6p (2.19) 

where {q){= q) and {p){= p) are the quantum-mechanical averages and 5q, Sp are the 
operators. By construction {6q) and {6p) are zero and [Sq,6p] = ih. Using Eqs. (j2.18j) and 
()2.19|) in {V'{q)) and a Taylor series expansion around (g) it is possible to express Q{t) as 

^-^ nl 

n>2 
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Here V^"'\q) is the n-th derivative of the potential V{q). For example, the second order 
Q(t) is given by Q(t) = -lV'"{q){6f). The calculation of Q(t) 20i2i, 22,23,24,25,26,27,28 therefore 
rests on quantum correction terms, ((5g"(t))(= Bn(t)) which are determined by solving the 
quantum correction equations as discussed in the Sec. IV. 

The c-number Hamiltonian corresponding to Langevin equation ()2.15t I2.16|) is given by 



2/i 



1 1 ^ r 2 ^ 

vi,Hj:-v"'H,mt) ,j:{^yfi....^^,f] (2,21) 

n>2 J i=l 

Note that the above Hamiltonian is different from our starting Hamiltonian operator 
()2.1|) because of the c-number nature of 1)2.211) . {xi,pi} are the quantum mean values of the 
co-ordinate and the momentum operators of the bath oscillators. 

To characterize the properties of the bath we define, as usual, the spectral density function 

as 

N 2 

J{oo) = ^y^^S{u;-u;,) (2.22) 
2/i ^—^ rriiUii 

Splitting the potential into a linear and nonlinear part as 

V{q) = ^fxuy + V,{q) (2.23) 

where Vi(g) is the nonlinear part of the classical potential V{q) we express, using ()2.23p . 
the quantum correction term as 

^i^V-W(g)5„(t) = l^^B,{t)uj', + V2{q) (2.24) 

n>2 

n>2 

Therefore the c-number Hamiltonian (Eq. 2.21) can be written as 

H = Ho + VM (2.25) 
where, Hq, the harmonic part of the Hamiltonian is given by 
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and VAr(g) is given by 



VN{q) = VM + V2{q) (2.27) 

V2{q) is the quantum correction term due to nonlinear part of the system potential. We 
define the mass weighted co-ordinates as 



"3 V'^i -^3 



q = y/JI q and x'j 
Diagonalizing the force constant matrix T of the Hamiltonian (Eq.2.26) as 



(2.28) 



UT = \^U (2.29) 
where U provides the transformation from old co-ordinates to the normal co-ordinatest^^ 



P 

yi 
y2 



( a' \ 



U 



X 1 



X2 



(2.30) 



\yN / \x'n J 

Here p and r/i are normal co-ordinates corresponding to system and bath modes respec- 
tively. 

The normal mode Hamiltonian for the harmonic part Hq is then given by^, 



N 



H 



NM 



(2.31) 



2' 2 

The eigenvalues Aq and are expressible in terms of the coupling constants of the system 
and the bath implicitly as follows: 



A^ = ujII 
A? = U31I 



1 ^ 



i = 1,2. ..AT 



(2.32) 
(2.33) 
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where fl2.32|) and ()2.33|) correspond to normal mode frequencies of the system and i-th 
bath oscillator respectively. 

The transformation (|2.3Up explicitly implies, 

N 

q' = Uoo P + ^ Ujo Vj (2.34) 

i=i 

where the matrix elements Uqq and Ujq can be expressed as 



u, 



00 



u. 



iO 



1/ 
1/ 



1 ^ 



1,2.. .AT 



(2.35) 
(2.36) 



For the present problem of dephasing it is significant to consider the coupling between 
the system and the bath modes to be weak. One can make use of the perturbation technique 
and Aq and A| are then expressible in simple terms aa^, 



Uij 



~t "^J ^1 (^1 - ^o) 



0(c 



1 + 



II rrij u;] {u] - uj^) 



+ 0(ct), J = 1,2 N 



l + 0(cj 



+ 0(c2), J = 1,2 N 



(2.37) 



+ O(c2), ij^O 



III. QUANTUM VIBRATIONAL DEPHASING RATE FOR A CUBIC OSCILLA- 
TOR 

It has already been established^'^'^ that harmonic oscillator model is not sufficient for 
a quantitative description of vibrational dephasing rate. The essential paradigm for the 
theory that has been used over the decades involves cubic nonlinearity so that the potential 
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assumes the following form, 



(3.1) 



Here / is a constant due to cubic nonlinearity. With Eg. ()3.1|1 the full Hamiltonian 
iJ[ (l2.25p ] in normal co-ordinate is given by 



N N 

H = Hnm + kuiiuoo P + ^^Ujo yjY + 3fiB2ulQkui{uoo p + ^Wjo Vj) + P^^'^Bsul^km (3.2) 

i=i i=i 

where we have used Oxtoby's notation, km = f /Qjj?/"^ and relations ()2.27|) and ()2.28|) . 
Here the first term denotes the normal mode Hamiltonian for the harmonic potential and 
the second term corresponds to classical nonlinear part of the potential. In addition to 
a constant shift of quantum origin third term signifies the quantum corrections to system 
normal mode where nonlinearity and quantum effects are entangled. In what follows we 
show that this term provides a substantial contribution to the vibrational dephasing rate. 

The anharmonicity in the potential shifts the minimum and the frequency of the system 
normal mode so that by applying the usual condition 







(3.3) 



to obtain the instantaneous minimum of the potential, pe, we have 



Pe 



000 



/2 



N N 

-(Ag + 6^111^00 ^ UjQVj) + (Ap + 12Ao/ciiiMoo ^ Uj^Vj - 36/i52Moo^iii)^ 

(3.4) 



The instantaneous frequency is given by 



Ao(t) 



1/2 



Pe 



An 



An 



1 + 



12 Unn ^111 



N 



^00 
'^O 



36 p B^it) kl 



11 



i=i 



\4 



1/4 



3 fciii 9 p B2{t) 

K 





(3.5) 
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where we have used uqo in the leading order. The instantaneous frequency shift is there- 
fore, 



Auj{t) = Ao(t) - Ao 



N 



3kni ^ 9/iS2(t)fci\i 



A, 



^0 



(3.6) 



In the weak couphng hmit the dephasing rate is expressed as, 



ep 



(Acj(t)Acj(O)) dt 



(3.7) 



where the averaging is carried out over the thermally distributed bath modes 



(A^(t)A.(0)) = %j:^ 



^0 -I L 



-nXj coth , ^ 



cos(Ajt) + 



11 



S2(t)52(0) (3.8) 



Here we have used the relations for the thermalized bath mode a^^i^^i^-^i^^i^^i^^i^^i^^ . 



{y^ms = {mms = {y^m^ms = 



2kBT 



(3.9) 



The quantum dephasing rate is given by 



9kf 



11 



X2 
-^0 



poo N 2 

/ 

•'» i-i 



-/lAo coth 



, 81/i2A;fii52(0) , , „ , , 
cos(Ajt) + ^ ^ ^ / B2{t) 

An 



(3.10) 

Using ()2.37|) for the expressions Aq, A^ and in ()3.10|) we obtain in the weak coupling 



regime 



^dep 



9kl 



2 /-oo 
11 



^0 



^ /i ruj J] [J] - ^2)2 



1 / hhj^ 

-nw^ coth , ' 

2 ^ V2A;bT 



cos(u;jt) 



Use of Eq. ()2.22|) in the above expression and continuum limit results in 



(3.11) 



11 



is^iii r, r, Ji^) 1 



Udep = T- / dt du 



<^P ~, ,2 / / ""^ , , (, ,2 , ,2\2 



1 , / huj 

—TujJ cotn 



2 V2A;i3T 



cos(u;t) 



oo 



+ ^ ^6 / ^^^2(t) (3.12) 

^0 ^0 

This is the general expression for quantum vibrational dephasing rate. The essential 
content of this formulae in addition to the usual first term obtained earlier by Levine et al^ 
is the second term signifying the quantum contribution to dephasing rate arising out of the 
nonlinearity of the system potential. This term is independent of the quantum nature of 
the thermal bath. An evaluation of this term requires the explicit calculation of the integral 
over quantum correction term B2{t) which we pursue in the next section. Keeping in view 
of the fact that J{uj) does not involve a specific choice of form for density of bath modes, 
we find that the expression for the dephasing rate as derived above is fairly general. 

The above method is based on the normal mode Hamiltonian of Pollak^^ adopted to 
a c-number description. An analysis of pure dephasing of a nonlinear vibrational mode 
has been worked out earlier to calculate non-Markovian line shape by Georgievskii and 
Stuchebrukhov^ using normal mode Hamiltonian treated by thermodjTiamic Green's func- 
tion approach. While the basis of present calculation of dephasing rate is Eq. ()3.7|) . the 
authors of Ref.5 have taken recourse to a different strategy to calculate the line shape. The 
differences in formulations and starting Hamiltonians (In Ref[5] a Leggett-Caldeira form of 
Hamiltonian, i.e., Eq.(j2^ without a counter term has also been employed) notwithstand- 
ing, the effect of a quantum contribution to dephasing width related to anharmonicity of the 
oscillator has been calculated in both Ref[5] and present analysis. The effect is due to the 
fact that the frequency of fundamental transition of a quantum nonlinear oscillator differs 
from harmonic frequency. To this end a continuation of the present analysis to calculate the 
associated frequency shift is instructive for comparison with those of others^. For this we 
return to the expression ()3.5|1 for instantaneous frequency Ao(t), which after keeping terms 
upto kill may be written as 



Ao(t) = Ao 



1 + 



Sfciii^ 21kl,J^ y 9/xi?2(t)fcLi 
2^ «io vAt) - I 2^ UjQ y.j{t) 1 

j^i \j=l / 



(3.13) 



where we have put Uqq to the leading order (unity). The time average frequency Aq is 
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given bj*^ 



An 



lim — 

r^oo T 



Ao(t') dt' 



Putting (prn^ in (pm|l we obtain 



(3.14) 



Ao — Ao 



^Aj 
2kBT 



9 fi B2 A^iii 



^0 



where we have used fl3.9p and B2 (= Uc Jq^"^" B2(t')dt' , being the cutoff frequency) 
is given by (the exphcit form of B2{t) is calculated in the next section). Furthermore 
with the replacement of Ujq using Eq. ()2.37|) and then using Eq. ()2.22j) in the continuum limit 
we obtain 



(3.15) 



Ao = Aq - ^ ^"^^ / J{uj) COth I :r^^ 1 duj ' 



2 vrAot^o Jo 



2kBT 



Since Ag is given by^, 



Ao = ^0 + - 



TT 



J(^) 



a; 



one may use ()4.9p to obtain 



(3.16) 



(3.17) 



Aq — 1^0 



1 + 



7^c 



The frequency shift Auq is then derived from Eq. ()3.16|) and Eq. ()3.18p as 



(3.18) 



AUJQ = Xq-UJo 

= ^ - r fivcoth f ^) - HMn^ (3,19) 

It is important to note that while in addition to the first two terms corresponding to 
treatment of Levine, Shapiro and PoUal^ the last one refers to temperature independent 
anharmonic quantum contribution proportional to kl^^ responsible to the frequency shift 
noted earlier in Ref . [5] . We mention in passing that the presence and absence of the counter 
term in the Hamiltonian may cause a significant difference in frequency shift with respect 
to direction towards blue or red region^. 
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IV. CALCULATION OF THE QUANTUM CORRECTION DUE TO NONLIN- 
EARITY OF THE SYSTEM POTENTIAL 



It has already been pointed out that a leading order quantum correction due to nonlin- 
earity of the potential of the system provides an important contribution over and above the 
usual expression for dephasing rate. To calculate this term explicitly we now return to the 
operator equation ()2.2p and use ()2.18p and ()2.19|) to obtain 

fiSq = 5p (4.1) 



6p + f 7(t - t') 6p{t') dt' + V"{q) Sq + J2^- V^"+'\q) {6^ - (5^)) = F{t) - f{t) (4.2) 

■^'^ n>2 

We then perform a quantum mechanical averaging over bath states with Ililillo^ilO))} 
to get rid of the term F{t) - f{t). The Eqs.dH)) and along with (jTTH|l and ^TH^ 

form the key element for calculation of the quantum mechanical correction. Considering 
the friction kernel 7(t) to be arbitrary (but decaying) we may calculate the leading order 
quantum correction for the harmonic mode for which higher derivatives of V{q) in (|4.2|) 
vanish. Now Eq. ()4.2|) becomes 



5p(t) = - -lit- t') 5p(t') dt' - fi ujl5q{t) 



(4.3) 

where fi uJq = V"{q) corresponding to the harmonic mode. The above equations ()4.1|) 
and (|4.3|) can then be solved by Laplace transformation technique to obtain 



6p{t) = - 6p{0) C,{t) + 6m C,{t) 



(4.4) 



where 



C,{t) = L 



-1 



+ s 7(s) +uJq 



(4.5) 



and 



C,{t) = l-uj'o / C,{t')dt' 
'o 



(4.6) 
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and 7(5) is the Laplace transform of 7(t) defined as 7(5) = l(t)e ^^dt. After squaring 
and quantum mechanical averaging Eq. ()4.4|) yields 



{Sfit)) = \ {5p\0) )Cl{t) + (5g2(0) )Cl{t) + - C^{t) C,{t){5p{0) 5g(0) + 5g(0) 5p(0)) 



For a minimum uncertainty state we chose^i 



(4.7) 



('^r(0)) = ^, W(0)) = ^ ((5p(0)5g(0) + (5g(0)5p(0)) = (4.8) 

Furthermore we assume the form of the spectral density function, J{uj), as, 

J(^) = \iu; (4.9) 
where 7 is the static dissipation constant. 

Using Eqs. ()2.4|l . ()2.22j) and (j4.9j] in the continuum limit we have 

7(t) = ^75(t) (4.10) 

Laplace transform results in 

7(s)=7i (4.11) 

where 71 = 7/2. 

Now with the form of 7(5) as given by Eq. ()4.1Hl the relaxation functions C„(t) and Cq{t) 
become 



1 

2cj^ 
2cJi 



S2 



-S2t 



(4.12) 
(4.13) 



where. 



nl/2 



71 



71 



+ CUl 



(4.14) 
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Making use of Eq. ()4.8|) . ()4.12p and ()4.13p in ()4.8p we obtain quantum correction term 
B,{t){= {Sem as 



1 + 



-2sii 



+ 1 + 



UJn 



-2s2t 



-4e 



(4.15) 



The above term can be utilized in the integral of the second term in Eq. ()3.12|) for its 
explicit evaluation to find out the dependence of the system parameters on the dephasing 
rate analytically. For better accuracy the systematic corrections to higher order can be 
worked out as discussed in detail in Refs. [18-24]. 



V. VIBRATIONAL DEPHASING RATE; COMPARISON WITH EXPERIMENTS 
AND DISCUSSION 

Having obtained the explicit expression for the leading order contribution B2{t) from 
Eq J4.15I we are now in a position to write down the total quantum vibrational dephasing 
rate. To this end we make use of Eq l4.15l in the second term and Eq. ()4.9|) in first term of 
the expression (j3.12p and obtain, after some algebra, 



with 



f^dep = Ki + K2 (5.1) 



= dt du — ^ coth — - cos^t 5.2 



and 



81 kf^^-f 
16 u 

The vibrational dephasing time can be defined as 



= (5-3) 



1 



(5.4) 

^dep 

The expression 1)5.11) is the central result of this paper. We already pointed out that K2 is 
a new contribution of quantum origin due to nonlinearity of the system potential, ki in the 
limit ksT ^ hwQ is the standard well known expression for the classical dephasing rate. It 
is important to note that Hi incorporates quantum effect due to heat bath only. Although 
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both Ki and K2 are dependent on nonlinearity, K2 vanishes in the classical limit and is also 
independent of temperature. The temperature dependence of the dephasing rate is due to 
the first term ki of Eq. (j5.2j) . It is important to note that at very low temperature as the 
integrand i^cj coth(/ia;/2A;sT) in Ki reduces to ^hu, the temperature independent vacuum 
limit, one observes that the dephasing caused by anharmonicity of the vibrational mode 
does not vanish even at absolute zero because of the contributions of these two terms. This 
aspect of temperature independence of the width of the transition from the ground state 
had been noted earlier in Ref . [5] . The origin of the temperature independence in K2 and the 
nature of dephasing may be traced to the second term in Eg. ()3.8|) which results from the 
third term of the normal mode Hamiltonian Eq. (|3.2j) . This term contains the quantum con- 
tribution to nonlinear potential explicitly calculated in terms of B2{t) in Sec. IV. This lowest 
order quantum fluctuation (or uncertainty) is independent of the quantum character of the 
heat bath and also temperature and causes frequency fluctuation leading to dephasing and 
therefore a homogeneous broadening of the transition similar to natural linewidth. In order 
to assess the relative contribution of the two terms in the total dephasing rate we estimate 
the numerical magnitude of these two quantities as well as the dephasing time for three 
diatomic molecules, N2, O2 and CO and compare them with experimental results obtained 
from either picosecond pump-probe technique^ or from Raman linewidth measurement of 
liquids using interferometric technique o^^i^^i^^i^^i^^ . We have also studied the rate as a func- 
tion of temperature away from critical point or triple point of these liquids and compared 
with experiments. 

The values of the parameters essential for calculation of the dephasing rate using formulae 
fl5.1|) are given in the Table 1. Apart from mass /i, frequency of the fundamental Uq, size r, 
at a temperature T two sensitive parameters are the static friction due to the liquid, 7, and 
the anharmonic force constant km. Although use of local viscosity which formally takes 
into account of wave vector dependence of the viscosity for the size of the probe has been 
advocated for calculation of friction, we confine ourselves to standard Stoke's expression (7 = 
6iTr]r/fi, f] being viscosity coefficient of the liquid). For diatomics we determine km from the 
spectroscopic constants and using^"^ km = — ^^3/^^ ^3 (1 + ^^)- The integrals in ()5.2p 
are calculated numerically for all the three cases. The magnitudes of ki and k,2 are shown 
separately in the table along with the percentage contribution of K2 in the total dephasing 
rate n^ep- Three pertinent points are to be noted. First, it is well known that classical 
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dephasing rate {i.e. ki in the classical limit) is higher than the corresponding quantum 
rate ki. This is somewhat reminiscent of the quantum supression of classical values of rate 
co-efficient for the thermally activated processes for shallow cubic potentials-^. Second, it 
is evident that just as in the classical theory, anharmonicity contributes significantly to 
the total quantum dephasing rate. In case of N2 it is as large as ~ 26.5 percent of total 
dephasing rate. The dephasing time thus calculated corresponds fairly to that obtained 
from experiments. Third, the quantum effects due to heat bath appears quite significantly 
through the integrals in ki. This is because, the frequency dependence of the integrand 
is quite sharp at around ujq as a result of the the frequency denominator (cj^ — ujq) and 
h ujq ^ ks T ioT these molecules. Therefore the diatomic oscillator behaves more closely as 
a quantum oscillator rather than its classical counterpart. 

The temperature dependence of the dephasing rate according to Eq. ()5.1|) has been com- 
pared for the set of parameter values mentioned in the Table 1 to that obtained from 
experiment o^^ 1 in Fig.l for the molecules N2, O2 and CO. While the experimental results 
span a wide range of temperatures covering triple point and critical point, we confine our- 
selves in the region away from the points of phase transition, since the present theory is 
outside the scope of phase transition. We find a linear dependence of dephasing rate on 
temperature which is a fair agreement with experimental observation a^^i^'^ . 



VI. CONCLUSION 



Based on a quantum Langevin equation and the corresponding Hamiltonian formulation 
within a c-number description, we have calculated quantum vibrational dephasing rate for 
a cubic oscillator system using first order perturbation technique and compared the results 
with experiments. It is shown that the vibrational dephasing rate comprises of two terms. 
The first term is the standard expression obtained earlier by several workers using inde- 
pendent approaches and is reduces to its classical limit as ks T ^ h ujq. This term is 
responsible for finite temperature dependence of the dephasing rate. The second term, an 
essential content and offshoot of the present theory, is a leading quantum correction due to 
nonlinearity of the system potential which contributes quite effectively to the total dephas- 
ing rate. Keeping in view of this observation, we believe that this quantum correction term 
is likely to play a significant role in triatomic and polyatomic molecules as well. 
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Figure Captions 



Fig.l: Variation of vibrational dephasing rate (FWHM) witli temperature (T) for the set 
of parameter values mentioned in the Table 1 and comparison with experiments for three 
different diatomics N2 (A experiments^ ; dotted line, theory), O2 (• experiments^; dashed 
line, theory) and CO (■ experiments^; bold line, theory). 
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Table 1. Quantum vibrational dephasing rate: comparison between theory and experiment for dia 

cOo T r li r\ km Ki Kj K^p 

Molecule (cm-^) ("K) (10"' cm) (10"'' g) (cP) (10^' g^'W^s"') (lo' s"') (lO' s"') (lO's' 

N2 2326 77 1.107 1.16304 0.158 7.47955 5.2965 1.90446 7.200 

O2 1552 90 1.345 1.32824 0.190 3.01 14.5996 3.6535 18.253 

CO 2170 77 1.214 1.13904 0.17157 8.1601 11.6269 6.5691 18.19 



^Present theory 

''Oxtoby's classical theory [Ref.(3)and (4)] 
"^Experiment [reference (12) and (14)] 
'^Experiment [reference (15)] 



